Spin torque contribution to the ac spin Hall conductivity 
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Using the recently proposed definition of a conserved spin-current operator [J. Shi et al, Phys. 
Rev. Lett. 96, 076604 (2006)] we explore the frequency dependent spin Hall conductivity for a 
two-dimensional electron gas with Rashba and Dresselhaus spin-orbit interaction in response to an 
oscillating electric field. We show that the optical spectrum of the spin Hall conductivity exhibit 
remarkable changes when the new definition of spin current is applied. Such behavior is mainly due 
to a significant contribution of the spin torque term which is absent in the conventional form of 
the spin current. In addition, it is observed that the magnitude and direction of the dynamic spin 
Hall current strongly depends on the electric field frequency as with the interplay of the spin-orbit 
coupling strengths. 



I. INTRODUCTION 

The spin-Hall effect (SHE) is a phenomenon that has 
and still motivating a very conspicuous interest among 
the semiconductor spintronicsi^ researchers. The SHE 
refers to a spin-accumulation induced via a driven spin 
current in response to a perpendicular dc electric field in 
finite width nonmagnetic systems which experiment spin- 
orbit interaction (SOI) effects^^ This spin (Hall) accu- 
mulation with opposite spin-magnetization at the edges 
of such systems has been beautifully demonstrated by re- 
cent optical experiments.-'- These observations were fol- 
lowed by purely electrical measurements of the SHE in 
metallic conductors^ 8 - and by its very recent detection at 
room temperature via Kerr spectroscopy in n-type ZnSe 
samples £ 

The role of disorder in the universal characteristic 10 
of the dc spin-Hall conductivity (e/8iv) in two dimen- 
sional electron gases (2DEG's) with Rashba SOI has been 
a central topic. Nowadays there is a general consensus 
that the spin-Hall conductivity is suppressed only in the 
static limit even in arbitrary weak disordered systems 
and as long as fc-linear Rashba SOI is considered in in- 
finite size samples . 10 i 11 ' 12 i 13 ' 14 i 15 This however does not 
hold in the presence of magnetic fields and/or magnetic 
impurities ) 10 ' 16 ! 17 neither for 2D hole system s 18 ' 19 ' 20 for 
which the static spin-Hall conductivity is robust against 
disorder. Recently Grimaldi et al£^ have showed that 
this is also true when the Fermi energy is comparable to 
the spin-orbit splitting energy (Ep ~ A). 

Yet another possibility to gain further physical insight 
to the SHE is the finite frequency regime of transport. 22 
Recent studies of the ac field-induced charge 23 and spin 24 
response of 2DEG's have emphasized the importance 
of the dynamic regime^ It has been suggested that 
an ac probing field can be used to control the spin- 
Hall current in 2DEG's with Rashba and/or Dresselhaus 
SO L 24 ' 25 It turns out that for finite frequencies (in the 
THz range), the cancellation of the intrinsic spin Hall 
effect due to impurity scattering 26 is no longer perfect, 
and in principle the effect may survived 

Sugimoto et al. 2S , on the other hand, have raised the 



question recently about the dependence of the vanishing 
of the spin-Hall conductivity cr s #(0) on the actual defini- 
tion of the spin-current operator used. This is certainly 
a fundamental issue that we extend here at the finite 
frequency regime of the spin-Hall conductivity, <j s h(uj). 
The widely ad hoc form used for the spin current op- 
erator for an electronic system J s = j{a z ,v y }, where 
a z is the Pauli spin z— component and v y the electron 
velocity operator, respectively, has the desirable form 
that resembles the charge-current operator. In 2D holes 
gases the conventional definition seems to be consistent 
with edge-spin accumulation experiments^ 9 - as with mea- 
surements of optically injected spin-polarized currents 
in semiconductors. 30 However, it does not fulfill a sim- 
ple continuity equatio n 31 ' 32 ' 33 and appears as incomplete 
in describing spin transport in systems with spin-orbit 
coupling^ 

Many efforts to clarify this issue have gone into this di- 
rection MelyM^ML^^^LB It seems nevertheless 
an unsettled problem yet. Furthermore, the discussions 
have been focused namely in the zero frequency limit. 
In particular, the recent proposal by J. Shi et al£^ of 
a dubbed unambiguous and proper definition of a "con- 
served" spin current density, is indeed physically appeal- 
ing and certainly deserves a close examination. Such new 
definition adds to the conventional part 3 S , a spin source 
term (torque dipole density P r ) associated to the elec- 
tron spin precessional motion with the total spin cur- 
rent density defined as J s = J S +P T . Moreover, among 
other interesting properties, it ensures that a continu- 
of- + V • Js = is always satisfied, with 



ity equation + V 



S z = ^^a z ^ describing the spin density. It is thus 
worth to elucidate to what extent the use of this con- 
served spin current operator provides new insight on the 
ac spin-Hall conductivity response in 2D systems with 
SOI effects. Whether or not such new definition of the 
spin current operator offers a physically satisfactory de- 
scription of the ac spin Hall conductivity should be ulti- 
mately validated by the experiment. 

In this paper we study the ac spin-Hall conductivity 
within the linear response theory employing the new def- 
inition for the spin current operator reported by J. Shi 
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et ali^i for a 2DEG in the presence of Rashba and Dres- 
selhaus SOI and subject to an alternating electric field of 
frequency to. It is verified post facto that the optical spec- 
trum of the spin Hall conductivity changes drastically 
when the new definition is used due to a non-trivial spin- 
torque contribution and to the interplay of the Rashba 
and Dresselhaus SOI. We predict that the torque dipole 
contribution to the ac spin-Hall conductivity generally 
dominates over the conventional part in typical samples. 

The work is organized as follows. In Sec. II we de- 
scribe the Hamiltonian model for a 2DEG in the presence 
of Rashba and Dresselhaus spin-orbit interaction. The 
main features of the Kubo formula in linear response ap- 
plied to the spin-Hall conductivity are presented in Sec. 
III. An analysis of the calculated ac spin-Hall response 
using the conserved spin-current operator is presented 
also in this section. In Sec. IV we give a discussion of 
the analytical and numerical results obtained for typical 
2DEG systems. A Summary is given in Sec. V and we 
end with a brief Appendix which outline the derivation of 
the spin-current-charge-current correlation function used 
in the Kubo formula. 




FIG. 1: (Color online) (a) Schematic diagram showing the 
angular anisotropy of the spin-split conduction bands at the 
Fermi energy in a 2DEG with finite Rashba and Dresselhaus 
spin-orbit interaction. The dotted C r (u>) curve is a rotated 
ellipse resulting from the condition e+(k) — £-(k) = hu>, see 
text for details, (b) Diagram depicting the allowed optical 
transitions (shadowed region) between the initial e_ and final 
e+ spin-split branches at a fixed polar angle 9. 



II. THE MODEL 



We consider a two-dimensional (2D) electron system 
with a single-particle Hamiltonian given by 



H 



2m* 



(1) 



where the spin-orbit part H so is the sum of the Rashba 
and (linear) Dresselhaus SOI type with 



a (3 

H SO = -T{0-yP X - a X Py) + —{0 X P X ~ CFyPy) (2) 

here, a and (3 are the coupling strengths for the Rashba 
and Dresselhaus couplings, respectively, with p x , p y the 
components of the 2D momentum operator and a x , a y 
the usual spin Pauli matrices. For the full Hamiltonian 
([T|) the energy spectrum is simply 



h 2 k 2 
2m* 



fikA(9) 



(3) 



with fi G {+1, — 1} denoting the spin-states for spin-split 
dispersion bands and A (9) = -J (a 2 + (3 2 )\(9) describes 
the angular anisotropy of the spectrum (see Fig. 1) with 



X(9) 



1 



2a[3 



■sin26l 



(4) 



and 9 = tan" 1 ( -jjjJM . When a or (3 is null, A(0) = I and 

the spin splitting in k— space becomes isotropic. The 
eigenstates have form 



fe(r)) = 



>2A 



1 

-/j,e~ 



(5) 



with the wave vector k= (k x ,k y ) = k(cos9,siri9) in po- 

— fiky 



lar coordinates, 6 — tan 



and A the area of 



the system. At zero temperature, the two spin-split sub- 
bands are filled up to the same (positive) Fermi energy 
level £f but with different Fermi wave vectors kp{9) = 

with 



+ k 2 \{9)-^k S0 ^{S)- Here, e F = -|^f 



kp = yj2im — 2k 2 Q and n being the electron density. The 
SOI introduces a characteristic SOI energy and wavevec- 

tor given by e so — r ,.__ s .° and k x ■ 

tively. 



2m* 



respec- 



III. SPIN HALL CONDUCTIVITY AT FINITE 
FREQUENCY 

Within linear response to a weak (spatially homoge- 
neous) electric field of frequency to oscillating in the 
plane of the 2DEG, the frequency dependent spin-Hall 
conductivity can be calculated using the standard Kubo 
formula— in the limit q — * 0, 



i{ui+iri)t 



a *y {uj) ~ hA^+Tr,) Jo 
x X;/( £ /.)T=o<^,,.WI[^f'"(*),««(0)]|^(r))dt, (6) 
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where we have assumed zero temperature and non- 
interacting carriers. The parameter r\ > is here to 
regularize the integral, but it can be viewed (phenomeno- 
logically) as a measure of the electron momentum dissi- 
pation effects due to impurity scattering events or any 
other many-body interactions^ with an overall time re- 
laxation rate r = rj^ 1 and /(e^) is the Fermi distribu- 
tion function. Here we employ the conserved spin-current 
operator of Shi et al£^ written in the Heisenberg pic- 
ture, JT%' z (t). Unlike the conventional definition J;,., the 
new effective spin current J s is (i) in conjugation with 
a mechanical or thermodynamical force (i.e. Osanger's 
reprocity relations are well established), (ii) it has the 
desirable property that it vanishes for localized orbitals 
which anticipates a zero SHE for insulators, and (iii) 
predicts that the static value of spin Hall conductivity 
has an opposite sign to those reported using the con- 
ventional definition for Rashba and linear Dresselhaus 
SOI modelsj 34 ! 39 It is shown in Ref.[34j], after imposing 
a nule torque dipole density outside the bulk, that the 
effective spin current operator J s can be simply recasted 
as J s = § Z ■ ^ s a resun % it adds an extra term 
which takes into account the spin torque contribution. 34 
In ([6]) v y (0) is the single-particle operator for the elec- 
tron velocity which can be obtained from the Heisenberg 
equation of motion v(0) = i [H, r]. For our system we get 

%(°) = -St-( ££d r £a )- The spin conductivity of/ (w) 
in ([6]) describes a z-polarized spin-current flowing in the 
x-direction as a function of frequency to in response to a 
weak electric field E(lo) chosen along the y-direction. 

To obtain explicitly J,^ z (t) we first take the x- 
component of the current operator with a spin moment 
polarized along the z direction at time t = (via the 
Heisenberg equation of motion) , yielding 



= ~~^ azPx+ {*> ~2h^ axPy + ayPx ^ 



-^(<?xPx + VyPy)} (7) 



2h 

where the symbol {,} denotes the anticommutator. The 
first term in Q represent the conventional spin-current 
part whiles the second term arises entirely from the spin 
torque contribution. 

It can be shown that the expected value of the 
spin-current-charge-current correlation function in (J5|) 

(^A r )\[Jx' Z (t),Vy( )]\^A r )) = J>(M) takeS the 

form (see Appendix) 



.7>(k,f) = ifi 



hkl (/3 2 



m*kA(9) 

2kA(9) 



t sin 



2kA(9) 



(8) 



By using this result in ([6]) together with the fact that 



J2 k e[{k-k+)(kp-k)], where 6(z) 



is the Heaviside step function, it is then possible to de- 
compose the frequency dependent spin-conductivity into 
the sum of two terms, namely 



in which the first one is given by 



(9) 



e(j3 2 ~a 2 ) f 2 * ,„ cos 2 « r*A°) 



de- 



47r 2 m* y A(0) J k+F{0) 

k 2 



dk 



[uj + i v ) 2 - {2kA(6)/h) 2 



(10) 



and comes from the conventional part of spin current 
definition. The second term is explicitly 



.H,.s_ e((i 2 -a 2 ) 



4ir 2 r 



271 in 8cos 2 9A(9) f k F {6) „ 
do — / dk 



[(uj + iri) 2 - (2kA(6)/K) 



212 



(11) 



and arises from the spin torque contribution to the net 
spin-current. As a general trend, it will be shown later 
that such contribution introduces a significant change of 
the spin-Hall optical response overcoming the conven- 
tional part for typical sample parameters. 

At this stage it is worthwhile to remark that in the 
limit lot — > oo the impurity scattering does not play a 
significant role in spin-transport. 10 Indeed, in the lan- 
guage of the diagrammatic technique, the resulting finite 
vertex corrections are negligible in the weak scattering 
limit for the ac high frequency field regime. This occurs 
because in linear response to the ac field, the perturba- 
tive expansion in powers of rii (impurity density) of the 
spin-current-charge-current correlation function goes as 
l/(wr) n which makes the contribution of the vertex cor- 
rections vanish at high frequencies i 10 ' 24 Thus the linear 
response Kubo calculation without including explicitly 
the vertex corrections should give a qualitatively good 
agreement with the full diagrammatic technique at high 
enough frequencies and low impurity densities^ 

The k— integration in expressions (fT0"|) and (TTTj) can be 
done analitycally, yielding for the conventional part of 
spin-conductivity 



1 a 1 



e/oiT tt a 
ftw 



cos 2 (9) 



1 



Se so X(0) 



tanh 



Se so \{6)h£a 



l&e F e so \{9) - H 2 Cj 2 



>d6 



(12) 



where Q = u) + irj and X(9) as is defined in eq. ((4]). The 
torque contribution becomes 
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T sH 



1 a 2 



e/8ir 
+ ■ 



7r a 



271 



tanh 



[16e F e so \{6) 



cos 2 f 
A(0) 1 
8e so A((9)fiu> 



^}d0 (13) 



[16e F e so X{9) - h 2 Co 2 } 2 - Ue 2 so \ 2 {e)h 2 Co 



Note that if j3 = (or a = 0), i.e. if A(0) = 1, the 
angular dependence of the integrand above reduces sig- 
nificantly and leads to a close analytic form for <j s c (ui) 
and a^. H (uj). For the case a ^ (3 the 8— integrals cannot 
be performed straightforwardly and a numerical integra- 
tion has to be implemented. Thus, it is convenient to 
consider first the physically reasonable limit k so <C kp, 
which typically holds for 2DEGs in III-V based semi- 
conductor heterostructures, and for which an analytical 
expression for the frequency dependent spin-conductivity 
can be derived. 

From (JTOJ) and (JTTJ) , we obtain to leading order in 

kso/kF> 



(0) 



1l/Vso - 1 
(V/Vso + I) 2 



(17) 



the Dyakonov-Perel spin- 



Since rj is 



where rj so = t so x , being t s , 
orbit relaxation time^ with r^ 1 = Qs^eJm 
typically smaller than rj so (i.e. loikp > hrj) note that the 
static value of a^ y z is always negative. 

It is known however that the dc limit is problematic 
within the Kubo formula when using the conventional 
spin-current operator, namely because it leads to the in- 
correct physics for e so hrj as a result of neglecting 
the contribution of the vertex corrections. This is not 
necessarily the case for finite frequencies and relatively 
low impurity densities (as discussed below), which is in 
fact the regime that we are mostly interested in here. 
In the opposite limit, e so Hrj, Kubo formula gives a 
result which coincides qualitatively with the expected re- 
sult [Hj]. We first consider this case. Taking the limit 
lo — > for finite but weak scattering mechanism with a 
time-relaxation rate t _1 = rj and expanding (| 15(1 and 
([To]) in powers of e so /hrj <C 1 to lowest order we get 



in which the conventional part reads 



(14) 



16(/3 2 - a 2 ) 



£f£so 



e/Sn (f3 2 + a 2 ) n^-^w 2 ] 1 ^' 

while the spin-torque part takes the form 



(15) 



o^M _ 32(/3 2 - a 2 ) 



e/8?r 



{H 2 



or 



£f£so 



rya-^ 2 ] 1/2 



eF£sofo 2 ti 2 [{6£F£so - h 2 oj 2 } 



ILK 



(16) 



with £ M = 2k F \ a + ji(3\. 

Our result JTSJ agrees with eq. (39) of Ref. [H 
obtained there via spin-susceptiblity calculation^ Note 
also that the second term of (fl"6|) vanishes for lo = 0, 
while the first one, when added to (fT5"]) . reverses the sign 
of the static value of the spin Hall conductivity (fT4|) . On 
the other hand, it will be seen that, at finite frequencies, 
the torque dipole contribution of Shi et al£^ generally 
dominates over the conventional spin current contribu- 
tion. 

It is illustrative to analyze the behavior of the spin- 
Hall conductivity in the static limit (lo = 0) in presence 
of weak disorder when the new conserved spin-current 
operator is applied. In particular, from eqs. (|15[) and 
([in]) , the static value, for the pure Rashba coupling case, 
can be written in the appealing form as 



gg(0) „, a 2 -P 2 [ 16£F£so 
e/87r ~ a 2 + I3 2 \ fi 2 rj 2 



(18) 



being cr^(0) = to first order. 

Now, if the impurity scattering is weak compared to 
spin-orbit coupling, e so /hri 3> 1, we obtain to first order, 



<(0) 
e/87r 



h 2 rj 2 (a 2 + P 2 ) 2 
me F £so (a 2 -(3 2 ) 2 



1 - 



sign(a 2 — (3 2 ) 



with a nonzero torque part 



(19) 



<(0) 
e/8vr 



ftV (a 2 +P 2 ) 2 



iepsso (a 2 -/? 2 ) 5 



sign(a 2 — (3 2 ) 



(20) 



The above expressions for cr^(0) reduce, in each case, 
to the known formulas for (3 — reported in Ref. [IF 



IV. RESULTS AND DISCUSSIONS 

We first study the case where only the Rashba cou- 
pling is present ((3 = 0). In Fig. 2 we plot the real part 
of the spin-Hall conductivity versus the frequency of the 
applied electric field as obtained from eqs. (fl~2|) and (fl3|) 
(A (8) = 1). The Rashba parameter has been fixed to a = 
1.6 x 10 -9 eVcm which is a typical experimental value for 
a 2DEG in InAs-based quantum wells. Using an electron 
effective mass of m* — 0.055 m , these values give us a 
characteristic Rashba energy of e so — 0.092 meV with a 
spin-splitting at the Fermi energy of An = 2akp — 5.6 
meV. The corresponding Fermi wave number is estimated 
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with kp = \Z2nn considering a moderated sheet electron 
density of n = 5 x 10 11 cm -2 . The parameter describ- 
ing the momentum relaxation rate has been chosen such 
that hi] — 0.4 meV, value that corresponds to high qual- 
ity samples with mobilities /i = er/m* ~ 5 m 2 /Vs and 
relaxation times of r ~ 1.6 ps. We have plotted as a ref- 
erence (dotted curve in Fig. 1) the result obtained using 
the conventional definition of spin-current [eq. (|10p]. 

The spin Hall conductivity a^ H (oj) shows resonance be- 
havior at energies hcu + ~ 2akp and foui- ~ 2akp, which 
correspond, respectively, to the minimum and maximum 
photon energy required to induce optical transitions be- 
tween the initial /i = — 1 and final /i = +1 spin-splitted 
branches (see. Fig.l). At low frequencies it approaches 
the universal value e/8ir, while it vanishes for high fre- 
quencies. Because of the finite value of the damping pa- 
rameter ?7, the spectrum has a smoothed shape across 

No- 



the shifted resonance frequencies w± = y lu± + rj^ 

tice that for k so <C kp the approximated solution of 
a^(uj) in Eq. (| 1 5 j) will give nearly the same numer- 
ical values with sligtly shifted resonance frequencies at 

where r\ > is un- 



;2±277V(Afl/fa7) s 



1. 



derstood and uj = rjyj (An/hrj) 2 + 1 is the intermediate 
frequency at which the conventional part vanishes. 

For the sake of contrast, the inset of Fig. 1 shows the 
behavior of the conventional part of the spin-Hall con- 
ductivity (eqjQ~Ql)) and that obtained from the expression 
(21) of Ref. [13J which incorporates explicitly the vertex 
corrections. It is shown that at finite high frequencies 
(hu > 0.5 meV) the result obtained neglecting the vertex 
corrections reproduce qualitatively the expected behav- 
ior, being in quite good agreement for finite frequencies 
outside the corresponding energy window referred above, 
that is for to > w^and lo < lu' + 

As for the effect of the torque dipole contribution on 
the spin Hall conductivity we notice that the first two 
terms in the expression (fT3|) change the sign of the terms 
giving a s c H {uj) in lfl2|). when added to obtain the total 
conductivity crf^ z (w). On the other hand, the last term 
of a^ H \ui) (ea. (|13p ) is a rational function which can be 
rewritten (for (3 = and (i^O) as 



~,2 



A|/n 2 )(a) 2 



4s|/ h 2 



{Co 2 - u\){Cj 2 - L0 2 _) 



(21) 



where er = 2e so = m*a 2 /h 2 . As is shown in Fig. 2, 
this contribution dominates the shape of the spectrum. 
It also resonates at the (shifted) frequencies uj±, being 
positive for u between uj' + and oj'_, and negative oth- 
erwise. Consequently, a dramatic change of the overall 
shape of the spectrum is observed. This is one of the 
main results of this paper. By measuring if possible, the 
spin Hall current /accumulation at low temperatures in 
the frequency domain, it could be helpful perhaps, to es- 
tablish the validity of the definition of the spin current 
operator proposed by J. Shi et al£^ by contrasting with 
our results. 
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FIG. 2: (Color online) Real part of the spin Hall conductivity 
as a function of the photon energy for a 2DEG system with 
Rashba SOI only (/3 — 0). The dotted curve is obtained us- 
ing the conventional definition of spin-current, whereas the 
dashed and solid curves are the torque contribution and the 
total spin conductivity, respectively, employing the conserved 
spin-current operator. In the inset the result obtained taking 
into account vertex corrections is contrasted with our calcu- 
lations. See text for parameter values used. 



New interesting features appear when the interplay of 
the Rashba and Dresselhaus SOI is considered. In Fig. 
3(a) we plot the frequency dependent spin-Hall conduc- 
tivity for P = 0.5a, here hq — 0.25 meV and the remain- 
ing parameters are as in Fig. 2. The results using the ap- 
proximated formulas (|15|) and (|16ll and the exact numeri- 
cal integration of (Tj"2")) and (TIB"]) are presented for compar- 
ison. Here is also evident the remarkable difference be- 
tween the optical spectrum resulting from the use of the 
standard and the new conserved spin-current operator 
definition, J^' z ■ In addition, the spectra become wider 
and highly asymmetric in comparison with those of the 
(3 = or a = case. As was recently discussed in Ref.l25l. 
the main spectral features can be understood as due to 
the anisotropic spin splitting caused by the simultaneous 
presence of the Rashba and Dresselhaus couplings. In the 
limit of vanishing temperature, the sum over states in eq. 
([5]) is restricted to the region between the Fermi contours 
<k <kp(0), for which e_(fc,0) < e F < e+(k 7 6), 
(see Fig.l). Certain distinctive frequencies can be identi- 
fied given the anisotropic fc— space available for the opti- 
cal response. To illustrate this, we consider the joint den- 
sity of states (JDOS) for the spin-split branches (Fig. 3b), 
which gives the number of direct transitions that can 
take place at the energy fuv. These transitions involve 
only states with wave vectors that satisfy the equation 
e+(k,6) — S-(k,9) — huo — 0, which for a fixed photon 
energy, defines a curve C r (u>) in A;— space describing a ro- 
tated ellipse with semi axis of lengths k a (ui) = hu>/2\a— j3\ 
and ki(u>) = Hu>/2\a+j3\ oriented along the principal axes 
(1,1) and (-1,1), respectively (see Fig.l). Thus, for our 
problem, the JDOS involves states only along the arcs 
of the resonance curve C r (uj) lying within the mentioned 
region kp(9) < k < kp(9). The peaks observed in the 
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FIG. 3: (Color online) (a) Frequency dependent spin Hall con- 
ductivity for a 2DEG system with finite Rashba and Dressel- 
haus SOI (/3 = 0.5a). Employing, the conventional definition 
of spin-current (dotted curve), the new conserved spin-current 
operator J^y Z in the limit k so -C k F (dashed curve) and by 
the exact numerical integration of eq. ((6]) (solid curve), re- 
spectively. Here hrj = 0.25 meV and other parameters as in 
Fig. 2. (b) Plot for the corresponding joint density of states 
(JDOS). There are four main frequencies here, two defining 
the optical absorption edges, lj_ and u+ , while the other two 
correspond to the peaks of the JDOS occuring at ui a and u>b- 
The latter two arises due to the symmetry of the spin-split 
conduction bands in k— space at the Fermi level. The dashed 
vertical lines at the frequencies referred above are here to 
guide eye. 



JDOS correspond to energy transitions involving states in 
the vicinity of the symmetry points k a (cu) = kp {n /A) and 
fcb(cj) = fc^(37r/4), for which the energy splitting reaches 
extreme values. These equations determine two energies 
tuu a = 2kp(n/4)\a-/3\ = 2k F \a-(3\+2m*(a-(3) 2 /ti 2 and 
hw h = 2fc+(37r/4)|a + /3| = 2k F \a + (3\ - 2m*(a + (3) 2 /h 2 . 
Similarly, we can see that there are absorption edges 
at energies fkJ± — 2k F \a =F (3\ T 2m*(a =F P) 2 /^ 2 , cor- 
responding to transitions between states at the points 
k a (ui) = kp(n/4) and = fc^(37r/4). For clarity 

such energies are indicated in Fig. 3 with dashed vertical 
lines. As expected, the spin Hall conductivity a^ z (tu) 
also shows structure at the photon energies hu>±, hu> a 
and hut- The finite value of fuq chosen here yields a 
slight shifting of these energies and an overall smoothing 
of the spectrum. 

It is clear that the new definition of spin-current opera- 
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FIG. 4: (Color online) (a) Spin Hall conductivity for f3/a = 
0.5,0.25,0.1 (solid, dashed and dotted curves, respectively) 
with the proposed conserved spin-current operator ' . The 
parameters are as in Fig. 2 with hr/ = 0.25 meV. Notice 
that the resonance peaks tend to separate in energy as the 
ratio /3/a is increased while its magnitude is diminished, (b) 
Color contour map of the ac spin-Hall conductivity showing 
its behavior with a continuous variation of the frequency u> 
and to the ratio (3/a. 



tor by Shi et al£^ yields dramatically different frequency 
response from that predicted by the conventional defi- 
nition. In addition, as it occurs with the pure Rashba 
(or Dresselhaus) SOI case, the torque dipole contribu- 
tion (not shown) turns out to be the dominant term in 
the spin-Hall conductivity. Notice that the spectra in 
Fig. 3 (a) show that the magnitude and the direction of 
the dynamic spin Hall current strongly depends on the 
frequency and on the spin-orbit coupling strengths a and 
/?, suggesting thus its control via electrical gating (by 
varying the Rashba coupling) and/or by adjusting the 
electric field frequency. 

We have also explored the effect induced of varying 
the ratio (3/a on the conserved spin-Hall conductivity as 
a function of the exciting frequency. In particular, in 
Fig. 4(a) the spin-Hall response is shown for the specific 
values of (3/a — 0.5,0.25 and 0.1 while fixing the rest 
of the parameters as in Fig. 3. Notice that the energy 
separation of the resonance peaks becomes larger and 
the intensity of the peaks get diminished as the aspect 



7 



ratio (3/a is increased. Such effect is emphasized in panel 
4(b) where a color map of the spin-Hall conductivity is 
plotted as a function of a continuous variation of the ratio 
[3 1 a and to. Such behavior responds to the widespreading 
{AS — ftio- — hw + ) of absorption bandwidth with 



V. SUMMARY 

In summary, we have examined the spin-Hall conduc- 
tivity in the frequency regime for a two dimensional elec- 
tron gas with Rashba and Dresselhaus spin-orbit interac- 
tion employing a recently proposed form for a conserved 
spin current operator. Our results shows that the opti- 
cal spectrum of the spin Hall conductivity changes sub- 
stantially when the new conserved spin current operator 
is used. It is predicted that the torque dipole contribu- 
tion typically overcomes the conventional part of the total 
spin Hall conductivity. In addition, it is shown that the 
magnitude and the direction of the dynamic spin Hall 
current is rather sensitive to the frequency and to the 
spin-orbit {a and j3) coupling strengths due to the spin 
torque contribution. We hope that these results will en- 
courage experimentalists to measure the spin Hall accu- 
mulation and/or spin density currents and to explore to 
what extent the new definition of the spin current opera- 
tor provides a satisfactory description of the ac spin Hall 
conductivity in such systems. 
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APPENDIX A 

Here we briefly outline the derivation of the expec- 
tation value for spin-current-charge-current correlation 
function (^k,p(r)|[j/'*(t),« y (0)]|^k, M (r)) = J>(k,t). 
We begin by writing the spin current operator at 
all times in the Heisenberg picture, J x ,z {t) = 
e lHt ' h J^ z {Q)e~ lHt / h , with j£' z {0) as given by 0. 



Similarly is done for the operators x(t), a x {t) and 
<J y (t), which after some algebraic manipulations explic- 
itly reads, 



V m h 



2A 2 p 2 

h_ 

2AV 



{(3 2 - a 2 )p y a z 
■W 2 -a 2 ) Py 



2Ap \ „ 
cos , _ i 

2tAp . / 2Ap 
sin 



h 2 



m 

(Al) 



K. 



<7 B (t) = <T X (0) + 



(ap x - (3p v ) 



K, 



Ap 

{aPx - (3p y ) 



AV 



2Ap , 
cos I -—7^t - 1 



rr- 



(A2) 



and 



<ry(t) = <J y (0) 



{(3p x - ap y ] 



■ ( 2A PA 



Ap 

(ftp x - ap y ) 



AV 



( 2A PA 



- 1 



(A3) 



with /C = a(a x p x + cr y p y ) — /3(a x p y + cr y p x ), and in which 
a and p are given in the Schrodinger picture. The ex- 
pressions above are needed in the calculation of the com- 
mutator [J x ' z (t),v y (0)] which gives 



_ (aa x + 0(r y ) iHt/h -wt/h 



2m* 



u z p x e 



Wt/h n „ r -iHt/n (aVx+Pvy) 
2m* 



2itp x 



< ''"OzPxe — - — — +i\x(0),C\-\ — -^C 

I > m* 



(A4) 



with C = -^{aa y - pa x ) sin(-^t) + -%P cos(^t 

and V — 2af3p x — (a 2 + (3 2 )p y - Calculating the expecta- 
tion value of (A4) using ([5]) we finally arrive to J-" M (k, t) 
as defined in equation (|SJ). 
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